arXiv:l504.05180vl [quant-ph] 21 Apr 2015 


Non—Contact Friction for Ion—Surface Interactions 


U. D. Jentschura^ and G. Lach^’^ 

^Department of Physics, Missouri University of Science and Technology, Rolla, Missouri 65409, USA 
^International Institute of Molecular and Cell Biology, Ksigcia Trojdena 4, 02-109 Warsaw, Poland 
^Faculty of Physics, University of Warsaw, Pasteura 5, 02-093 Warsaw, Poland 

Non-contact friction forces are exerted on physical systems through dissipative processes, when 
the two systems are not in physical contact with each other, or, in quantum mechanical terms, when 
the overlap of their wave functions is negligible. Non-contact friction is mediated by the exchange 
of virtual quanta, with the additional requirement that the scattering process needs to have an in¬ 
elastic component. For finite-temperature ion-surface interactions, the friction is essentially caused 
by Ohmic resistance due to the motion of the image charge moving in a dielectric material. A con¬ 
ceivable experiment is difficult because the friction force needs to be isolated from the interaction 
with the image charge, which significantly distorts the ion’s flight path. We propose an experimental 
setup which is designed to minimize the influence of the image charge interaction though a compen¬ 
sation mechanism, and evaluate the energy loss due to non-contact friction for helium ions (He"*") 
interacting with gold, vanadium, titanium and graphite surfaces. Interactions with the infinite series 
of mirror charges in the plates are summed in terms of the logarithmic derivatives of the Gamma 
function, and of the Hurwitz zeta function. 

PACS numbers: 31.30.jh, 12.20.Ds, 31.30.J-, 31.15.-p 


I. INTRODUCTION 

If an ion trajectory is parallel and close to a surface 
(close to a “half-space filled with a conducting material”), 
then it is intuitively obvious that the motion of the im¬ 
age charge inside the conductor will lead to some sort of 
Ohmic heating. Unless one can show that the motion of 
the image charge is frictionless, this Ohmic heating can 
only be compensated by a corresponding loss in the ki¬ 
netic energy of the projectile ion. Here, we devote special 
attention to a rederivation of the non-contact friction in 
ion-surface interactions, with an emphasis on a proposed 
experimental setup which serves to eliminate the “signal” 
(the quantum friction) from the “noise” (the interaction 
with the image charge). The latter can otherwise lead to 
a net gain or loss of kinetic energy for the ion. 

Provided one accepts the existence of non-contact fric¬ 
tion for ion-surface interactions, the same effect (but at a 
smaller scale) can be expected for atom-surface interac¬ 
tions. Namely, any atom constitutes, due to its quan¬ 
tum mechanical nature, an “oscillating dipole” whose 
frequency-dependent response to exciting radiation is de¬ 
scribed by the (dynamic dipole) polarizability. Because 
the dipole moment fluctuates, there is a fluctuating image 
charge moving in the medium, which also leads to Ohmic 
heating and by consequence, non-contact friction. How¬ 
ever, the friction force in this case will be much smaller, 
numerically, than in the case of ion-surface interactions. 

One can draw a distant analogy to ion-atom versus 
atom-atom interactions. It is well known that the ion- 
atom interaction (or more generally, the interaction of a 
charged particle with an atom) has the functional form 
1/R^ where R is the ion-atom distance. For two neutral 
atoms, the dominant term in the non-retardation limit 
is the van-der-Waals interaction, which gives rise to a 
l/i?®-interaction. For the quantum friction effect, with 


an ion-surface distance Z, the damping coefficient ry in 
the friction force F = —ryu is proportional to IjZ^ for 
ion-surface interactions, which is modified to a l/2^®-law 
for atom-surface interactions [1]. The pattern is clear: 
There are two more powers of the distance [R viz. Z) in 
the denominator for interactions with atoms, as opposed 
to the corresponding interaction involving an ion “at the 
other end”. 

Alternatively, we may treat the case of finite tempera¬ 
ture [2-6] as follows [1, 7-9]. When two physical systems 
are in contact, then the charge fluctuations in the one 
will induce mirror charge fluctuations in the other. At 
zero temperature, the alteration of the available modes 
for quantum fluctuations due to the presence of the sur¬ 
face (due to the “boundary conditions”) gives rise to 
atom-wall interactions of the Casimir-Polder type (see 
Chap. 8 of Ref. [10]). At finite temperature, one has to 
replace the integral over the virtual photon frequency by 
a summation over Matsubara frequencies [see Eq. (30) 
of Ref. [1]]. The zero-temperature limit is obtained by 
the replacement j3~^ Yl^=-oo ^da;/(27r), where 
the nth Matsubara frequency is = 27r/{h[3), and 
P = l/{kBT) is the Boltzmann factor. The same ob¬ 
servation is made in Sec. 81 of Ref. [11], where it is 
shown that in the zero-temperature and short-distance 
limit, the Casimir force between two solids reduces to 
an expression which is in agreement with the van-der- 
Waals force between atoms [limit of two dilute media, 
see Eqs. (81.1), (81.9) and (82.3) ff. of Ref. [11]]. The 
calculation of static forces is most easily accomplished in 
the imaginary-time formalism, which leads to the Mat¬ 
subara frequencies, while the calculation of friction forces 
is usually done using the real-time retarded fluctuation- 
dissipation theorem [1, 9]. Indeed, in Ref. [9], it is shown 
that the friction force on an atom in the vicinity of a di¬ 
electric surface can be obtained from a calculation of the 
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van-der-Waals interaction for an atom which undergoes 
a small periodic mechanical oscillation, after subtracting 
the conservative static, and the conservative oscillatory 
component (vibration component) of the van-der Waals 
force [see Eq. (24) of Ref. [9]]. 

In order to calculate the friction force, one has to de¬ 
scribe the correlation of electric field fluctuations at dif¬ 
ferent points, in the presence of a dielectric material fill¬ 
ing the half-space z < 0, due to charge fluctuations in the 
dielectric [2-4, 11]. We here use the real-time formulation 
of the fluctuation-dissipation theorem and apply it from 
first principles. Another point is that one has to apply the 
Wick theorem to the thermal fluctuations, remembering 
that the Wick theorem holds both for quantum as well as 
for thermal (statistical) fluctuations. Finally, the “ther¬ 
mal factors” need to be taken into account accurately, 
with a full account of the quantum nature of the prob¬ 
lem. Last, but certainly not least, one has to be careful 
in applying the conventions for the Fourier transform in 
the theory of thermal fluctuations and in electromagnetic 
signal theory [2-4] correctly. These can otherwise lead to 
inconsistent prefactors in the final results. 

We continue in Sec. II with a discussion of ba¬ 
sic concepts underlying the fluctuation-dissipation the¬ 
orem, which is central to the derivation of the finite- 
temperature non-contact friction effect, and with the 
thermal correlations of the electric field. The ion-surface 
interaction and the calculation of the non-contact fric¬ 
tion force on a charged particle are discussed in Sec. III. 
Finally, a sketch of a proposed experimental setup is dis¬ 
cussed in Sec. IV. It involves two parallel conducting 
plates which give rise to a series of mirror charges and 
requires the calculation of interaction potentials with the 
mirror charges. Conclusions are drawn in Sec. V. 


II. FLUCTUATION-DISSIPATION THEOREM 


A. Fluctuation—Dissipation Theorem 


We start by recalling the Boltzmann factor /3 and the 
bosonic thermal occupation number n(a;). 


P 


1 


n{uj) 


1 

exp(/3 huj) — 1 ' 


( 1 ) 


The Fallen-Welton thermal factor 0(w,T) is given by 
the relation 


Q{uj,T) = hoj (i -I- n{oj)) = ^hoj coth (i h/3uj) . (2) 

It has the properties, 

0(u;,T)= -0(-o.,T), (3a) 

(I-bn(u;)) (1-bn(-a;)) = ^ , (3b) 

(1 -b n{u})) n{uj) = - ^ ■ (3c) 


In the high-temperature limit, we have 0(a;, T) —>■ 1//3 = 
ks T. Let now x be an observable of a dynamical system 
with Hamiltonian Hq{x) subject to thermal fluctuations. 
We assume that x{t) fluctuates around its mean value 
(x)o with fluctuations characterized by a power spectrum 

S^{u}) = {x{t) a:(0))^ ^ ^ 

Let us consider the conjugate variable of x, namely, a 
scalar force field / which alters the Hamiltonian to 

H = Ho{x) + x{t) f{t). (5) 

The response of an observable x to the field term is char¬ 
acterized (to first order) by the susceptibility of the linear 
response function x(t) of the system, 

t 

{x{t)) = {x)o+ J x(t - t) /(r) dr , (6) 


where the mean value {x)q is obtained as a thermal av¬ 
erage taken over the distribution governed by the un¬ 
perturbed Hamiltonian Hq{x), i.e., with respect to the 
weight function 


Wo{x) 


exp [-j3HQ{x)\ 

J dx' exp [—/3iJo(a:')] 


( 7 ) 


The field term is adiabatically switched on at r = — c». 
Expanding the weight function of the full Hamiltonian 
given in Eq. (5) for small /(t), one can easily motivate [12] 
that the the imaginary part of the complex susceptibility 


X(a;) = Re [x(u;)] -b i Im [x(u;)] , (8) 

is related to the power spectrum of fluctuations, which 
by the fluctuation-dissipation theorem [2-6] reads as 

{x(t) x(O)),^ = Sa:{uj) = Iin , (9) 

(jJ 

This involves the Kallen-Welton thermal factor. In the 
high-temperature limit, one may replace Q{uj,T) —>■ 
1//3 = ksT. We should note that the conventions for 
the Fourier transform used in the current work follow 
those commonly used in physics, 

/(a;)= I dte-‘/(i), f{t) = J ^e--‘/(a.), 

( 10 ) 

and these are different from those in electrical engineering 
(see Refs. [2-4]). 


B. Correlation Function for the Electric Field 

Let us try to motivate a formula for the thermal cor¬ 
relation of the electric field in the vicinity of a dielectric 
material, based on the thermal charge fluctuations inside 
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the dielectric. We first observe that the scalar potential 
<i)(r) is given as 

<I)(r) = J d^r'p{7^), (11) 


where p is the charge density. In free space, the Green 
function is given as G{f,r') = (47reo|r — In the 

vicinity of a dielectric wall, one has to modify the Green 
function as follows (non-retardation limit, see Ref. [1]), 


G(r, f') —>■ G(a;, r, t^) 


1 

47reo|r' — f] 


e(uj) — I 1 

e(w) -I- 1 47reo|f^ — r + 2n±(f- n_L)| 


( 12 ) 


We here ignore a possible path difference of the emit¬ 
ted perturbation (at and the incoming wave (at f*). 
The normal vector n± is the outward normal pointing 
away from the half-space filled with the dielectric mate¬ 
rial. Indeed, we shall need the Green function in the limit 
this limit, the path difference of the reflected 
and the emitted wave vanishes. Gorrections to this result 
due to the velocity of the atom in the x direction are of 
order Ziuv^jc^ and are negligible on the order of interest 
for the current paper. Here, Z denotes the ion-surface 
distance. 

Our formula (12) is in agreement with Eq. (6) of 
Ref. [1], which gives the fluctuations of the scalar in¬ 
stead of the vector potential, due to charge fluctuations 
in the vicinity of a dielectric wall. The 1 / (27r) prefactor in 
Eq. (6) of Ref. [1] is due to the different Eourier transform 
conventions used therein, which follow the conventions 
commonly adopted in electrical engineering. The corre¬ 
lation function of the fluctuations of the scalar potential 
reads, in full agreement with the general paradigm set by 
the fluctuation-dissipation theorem given in Eq. (9), 


(4>(r) ci>(r'))^ = [G(a;; r, r^)] . (13) 

UJ 


Here, p is the “fluctuating force”, whereas $(r) takes the 
role of the fluctuating signal. The correlation function for 
the electric field follows by differentiation. 


(E;,(f)E,(r'))^ 


- (V,$(r)V'4>(r'))^ (14) 

- Im [V,V'G(c.; r, r')] , 

2Q{uj,T) / e(c^)-l \ 

w ve(w) + 1/ 


X V,,V' 


^ 47reo|f^ — 'c + 2n_L(r • nj_) 


The minus sign is explained because negatively correlated 
charge fluctuations at two points along a reference line 
generate fluctuations of the electric field which are pos¬ 
itively correlated along the line joining the two points. 
Here, €{uj) denotes the relative dielectric function (rel¬ 
ative permittivity), which has to be multiplied by the 


vacuum permittivity eo, if one would like to obtain the 
full electric permittivity of the medium. 

There is one last subtlety to discuss. We should be care¬ 
ful interpreting the correlation function 
Namely, according to Eqs. (32) and (49) of Ref. [I], it is 
more accurate to relate the following, symmetrized corre¬ 
lation function to the imaginary part of the susceptibility, 
according to the replacement 

/ CO 

dtcosiujt) {Ei{r,t)Ej{r,0)) 

-CO 

/ CO 

dt (e-‘ + e-'“‘) {E,ir,t)E,if^,0)) 

-CO 

/ CO 

dte“^‘ {E,{r,t)Ej{^,d)) 

-CO 

/ CO 

dte‘“* (E;,(r,-t)E;,(r“',0)) 

-oo 

/ CO 

dte“^‘ ({E;,(r,t),£'j(h',0)}) 

-OO 

= (i{E;,(r,t),E;,(r',0)})^ . (15) 

The distinction between {Ei{r)Ej{r')) ^ and 
(i{Ei(r, t),Ej(r',0)}) becomes important when 
more than two operators are involved [1, 13]. So, we 
should write, more correctly, 

{\{Ei{r,t),Ej{^ ,0)})^ = ^^^^Im 

* dTreol'f’’— r-f 2n_L(r • fi_L)| 

The result in Eq. (16) is consistent with Eq. (15) of 
Ref. [14]. 


e(^) - 1 

e(a;) -I- 1 


III. QUANTUM FRICTION FOR IONS 
A. Ion—Surface Interactions 


Armed with the results from Sec. H, we are now in the 
position to evaluate the friction force due to the mirror 
charge running in the dielectric, which generates Ohmic 
heating. The force on a charged particle is given by 

F{t) = ZeE{T,t), (17) 


where we assume e to denote the electron charge, and the 
ion is Z-io\d negatively charged. We also use the result 
[r = {x,y,Z) and f = {x\y',Z')] 


hnqVa:V(,— 

r'->-r \r 


1 

f+ 2n±{f- fi_L)| 


1 

8^ ■ 


(18) 


We consider the formula for the friction force. 


Fx = -rjVoc ■ 


(19) 
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According to Eq. (8.15) of Ref. [15] (for a more modern 
perspective Ref. [16]), the fluctuation-dissipation theo¬ 
rem determines the friction force via a Green-Kubo for¬ 
mula. In the derivation, we use the fact that for t > 0, 
we have {Fx{t) Fx{0)) = {Fx{0) Fx{—t)) due to time- 
translation invariance. We can thus symmetrize the inte¬ 
grand as follows. 


poo 

r] = l3 dt {Fx{t) Fx{0)) 

Jo 

= f dt{^{F4t),F,{0)}) 

= Jinr f dt E^{f',0)}) 

2 r^r> J_^ 


= ^{ZeY lim 

2 ^ ' i.j-s-0 


LO \e(aj)-|-l 


X lim V.rVl 

r—¥r‘ 


^ ^‘iTreol'r'— f + 2n±{f ■ n±) 


P{Ze) 


lim 


-Im 


47reo \e(u;)-|-l/J 8Z^ 


:(w) - 1 


(Ze)^ 


lim < — Im 


327reo <^^0 (a; \e(w) + 1/ J 


z(uj) - 1 


^ • (20) 


=L 


The end result is 


77 = 

L = 


(Ze^L 
327reQ Z3 ’ 

lim I — Im 

til—)-o 1 uj 


U(^) + iyi 


(21a) 

(21b) 


We note that the drag force F^ = —rj Vx is independent of 
the mass of the projectile particle; it only depends on its 
charge state. The result (21) depends on the asymptotic 
shape of the Boltzmann factor n{uj) = l/(exp(/3/ku) — 1), 
which goes as l/{Phuj), for small w. The final evaluation 
of the limit L proceeds via a consideration of the low- 
frequency limit of the dielectric function of the material 
and crucially depends on the lowest electronic resonance 
frequency of the material. 

A final word on retardation corrections is in order. In 
the treatment of Casimir-Polder interactions of neutral 
atoms, the parameter governing the retardation effects 
is [10, 17] 




LU Z 
c 


UJ 


2 2 
a Trie c 

h 


( 22 ) 


where we indicate a typical value of an atomic transi¬ 
tion frequency uj (the fine-structure constant is a). Re¬ 
tardation sets in when the phase of the atomic oscillation 
changes significantly over the time it takes light to travel 
to the surface and back, i.e., when ^ 1 and thus 


a 


137 oo, 


Oo = 


a me c 


(23) 


Retardation changes the leading 1/Z^ interaction for 
short distances to a IjZ^ term at long range [10, 17]. 

For the ion-surface interaction, the leading conserva¬ 
tive term, for any distance, is given by the 1/2^ attractive 
interaction with the mirror charge, and, for the dissipa¬ 
tive friction term, by the 1/Z^ form given in Eq. (21). 
The estimate of the retardation corrections then has to 
proceed differently; it is necessary to evaluate the retar¬ 
dation correction to the (in principle electrostatic) in¬ 
teraction of the ion with its mirror charge. One should 
compare (i) the time it takes light to travel to the surface 
and back to the ion moving alongside and parallel to the 
surface, to (ii) the time it would light to travel back to a 
static ion. The relative difference is a measure of the re¬ 
tardation correction to the electrostatic interaction. The 
expansion parameter in this case is easily found to be 

, (24) 

c 

where v is the ion’s velocity. Because these fly-by veloc¬ 
ities are much smaller than c in typical atomic-beam or 
ion-beam experiments, we can safely ignore the retarda¬ 
tion corrections. Further retardation effects due to the 
frequency of the exchange photon can be ignored; the re¬ 
sult given in Eq. (21b) is formulated in terms of a limit 
for small frequencies a; —0. 


B. Material and Friction Force 


A rather general functional form for the dielectric func¬ 
tion of a material is given as [17, 18] 


00 . 1 

, Z , , an — ibnUJ 

;(w) = 1 + > ^^- 

^ -i7„w 


(25) 


The expansion coefficient is an amplitude which can 
be written as 


CLn — UJn , bn — ACyi ^n ■ ( 26 ) 

Recently, an approach with manifestly complex expan¬ 
sion coefficients has led to an excellent fit of the dielectric 
function of silicon over a wide frequency range [18]. The 
static limit of the dielectric function then reads as 

OO 

e(w = 0) = 1 + ^ Ae„ . (27) 

n—0 

The functional behaviour of e(a;) for small w crucially 
depends on the value of the lowest resonance frequency 
uJn=o, where we reserve the subscript n = 0 for the case 
of an excitation at zero resonance frequency wq = 0 and 
otherwise set its multiplying coefficient to zero in the case 
of the absence of such a resonance. We ascertain that if 
a resonance with wq = 0 exists, then e(a;) — 1 = 0(uj~^) 
for small angular frequencies. This is the case for metals 
and other conductors, where, by definition, the energy 
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FIG. 1: Basic proposed setup for the measurement of ionic non-contact friction. An incoming ion beam is velocity- 
selected and enters an interaction region with a compensating slit configuration, involving two parallel bars of 
(preferentially single-crystal) material, whose induced non-contact friction is to be studied in the experiment. A 
further free beam path of a few meters in length serves to eliminate conceivable contributions to the ionic energy 
loss due to imperfect alignment of the interaction region. Finally, the energy of the emerging ion can be measured 
after it has passed through an aperture. 


required to excited an electron into the conduction band 
vanishes. By contrast, if oq = 0 and the lowest resonance 
frequency of the material is wi ^ 0, then e(a;) = e(0) -I- 
0(u!), so that Im[e(a;)] = 0(uj) for small w. In both cases, 
the limit L defined in Eq. (21b) is finite. 

Let us first consider the case of a conductor. Then, the 
dominant contribution for small oj is given by the term 
with n = 0 in Eq. (25), 


e(uj) « 1 - 


ao 


uj(cu -I- iyo) 


= 1 - 




Lu(uj -I- iyo) 


(28) 


where cUp is the plasma frequency of the free electron gas 
constituting the conduction band electrons. This leads to 


consideration of the lowest resonance. In the case bi = 0, 
we have 


e(uj) « 
L = 

1 = 


^ ' O 9 • ’ 

(Jjf — 

2ai7i 

(oi -b 2a;f)2 ’ 

{ZeY ai 7i 

16 TT (oi + 2w^)^ 


(32a) 

(32b) 

(32c) 


Here, we assume that the lowest resonance dominates, 
i.e., that ^ ujf 1^ oi. If a modified Lorentz profile 
with 7 ( is chosen (according to Ref. [18]), then the result 
simplifies to 


L = lim 

o;—>-0 



LU 


le(..) + i;/ 


270 


UJ„ 


2eo 

ctt(O) ’ 


(29) 


where 


ariO) = Co — 
70 


(30) 


is the temperature-dependent conductivity (measured in 
inverse H m) of the material, for a direct current (alter¬ 
nating current of zero frequency). We recover the result 
given in Eq. (35) of Ref. [1] for a conductor. 


jZer 

16 7 r(TT( 0 ) 


( 31 ) 


The only temperature-dependence remaining in this re¬ 
sult is via the implicit dependence of the “static” conduc¬ 
tivity (Tt( 0) on the temperature. In the low-temperature 
limit, many materials become superconducting, which 
implies that their conductivity crT(O) diverges and the 
friction coefficient rj vanishes. 

For an undoped semiconductor and other materials 
such as gases, the evaluation of the limit L proceeds via a 


e{uj) 


1 + E 


-i7n^) 
- i 7„ w 


T ^ 2 Acn (-fn - In) 

■^i(2 + E^Ae^)2a;r 

_ {Zef Acn (7n - 7^) 

^ 16^Z3 4^(2 + E„Ae«^)^wr 


(33a) 

(33b) 

(33c) 


where in the last line both sums extend over m,n = 
1,..., 00 . Let us match these formulas with the known 
expressions for a dilute gas (bulk material), whose di¬ 
electric function is given by 


e(a;) = 1 + ^ E q,(^) 

^ eo 

n 

_ J Ny _/nO_ 

eo V eIq - ir„ [hw) - {hujY ■ 


(34) 


Here, Ny = N/V is the volume density of atoms, and the 
Eno are the excitation energies from the ground state to 
the nth excited state. The (dipole) polarizability is a(w), 
and the dipole oscillator strength for a transition from the 
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reference state k to the virtual state n is fnk = fnk 
where 


47re^ 

/nO ~ 2 -fcTn i TT9 


(2£+l)2 

<EE 


inYyimif-i) 


(35) 


Here, £ = 1 denotes the dipole contribution to the polariz¬ 
ability, where 2^ is the multipole order [19]. We sum over 
the magnetic projections m of the spherical harmonic and 
over the magnetic projections mk of the excited state 4>k- 
The identification of Eqs. (33) and (34) then proceeds as 
follows. 


— T 7 Tn 

a 



Ae„ 


NyfnO 

<^0En0 ' 
(36) 


These quantities can directly be used in Eq. (33). 


IV. “SNIPER” SETUP CONFIGURATION 
A. Mirror Charges for Parallel Conducting Plates 

We consider an experimental setup as given in Fig. 1 
and strive to calculate the interaction potential of the 
ion with the mirror charges in both the upper, as well as 
the lower slab, as well as the induced friction force. For 
simplicity, we shall consider, in the following, two con¬ 
ducting parallel plates. The upper surface of the lower 
plate is in the xy plane (2 = 0 ), whereas the upper plate 
is at 2 = a. The ion’s 2 coordinate is denoted as Z. The 
following analysis is in part inspired by Refs. [20, 21]. The 
(positively charged) ion generates two image charges in 
the two conductors, which in turn, by mirroring them 
against the respective other conducting plate, generate 
an infinite series of mirror charges, which can be identi¬ 
fied as follows. Iterative calculation of the positions of the 
mirror charges leads to are two “upper” series of mirror 
charges, negative ones at position 2na — Z [the distance 
to the reference point Z is 2(na — Z)], and positive ones 
at positions 2na + Z (the distance is 2na). There are also 
two series of mirror charges in the lower slab, positive 
ones at positions —2na + Z (with a distance of 2na), and 
negative ones at position — 2 (n—l)a — Z, whose distance 
to the reference point is 2 ((n — l)a-|-Z). 

Let us consider the limiting case of the ion being close 
to the lower plate at 2 = 0, i.e., the limit Z ^ 0. The 
dominant interaction potential will be due to the clos¬ 
est mirror charge in the lower slab, which is located at 
2 = —Z. The distance of charge and mirror charge is 2Z, 
and the interaction potential is —(Ze)^/( 167 reoZ), where 
the prefactor takes care of the distance 2Z, as well as the 
fact that the electric field is zero inside the conductor (see 
the Complement to Chap. 11 of Ref. [22]). Alternatively, 
the additional factor 1/2 can be understood as follows: 


If one is to move the charge upward by a given distance, 
then the mirror charge moves “automatically”, and thus 
the work required to move the charge (and its interac¬ 
tion potential) is halved. The sum over all other mirror 
charges leads to the following total interaction potential. 


V = -^^^C{Z). 
47reo 


(37) 


A series representation of the correction factor is as fol¬ 
lows (again in the non-retardation limit, see Ref. [21]), 


- 00 

C{2) = 5 E 


n—1 

1 

' 4 ^ 


1 


1 


2a{n — 1) -t- 2Z 2an — 2Z 


an 




^ [ 1 - - 
a 


2 7_e 


, (38) 


where je = 0.57721... is the Euler-Mascheroni con¬ 
stant. The logarithmic derivative of the T function and 
its generalization 'I'(”)(x) read as follows, 

d d" 

^(a;) = = ^^( 2 ^) ■ (39) 

We anticipate that we shall need 4'*^"’=^)(a:) in the follow¬ 
ing derivations. For small Z, the asymptotic expansion 
reads as 


in agreement with the consideration sketched above. 

For the “sniper” configuration to be discussed below, 
the setup sketched in Fig. 1 suggests to assume an (al¬ 
most) symmetric configuration, with the atom in the 
middle between the two conducting slabs, i.e., one ex¬ 
pands about the point Z « a/2. For perfect symmetry, 
we have 


= = = (41) 

\ 2/ a a a 

If we assume perfect symmetry, then the two nearest mir¬ 
ror charges are at a distance a from the ion, which, to¬ 
gether with the correction factor 1/2 of Ref. [22], would 
suggest a value of 2 x l/(2a) = 1/a for the C coefficient. 
The contribution of the remaining mirror charges reduces 
this result to the value ln( 2 )/a. 

We now consider the contribution of the remaining 
mirror charges to the friction. To this end, we first re¬ 
call the result from Eq. (21), which is initially valid for a 
single surface. 


{ZefL 
~ 327reo Z3 


-t V(2) 


[ZefL 

327reo 


D[Z). 


(42) 


Here, the variable Z in the expression for 77 denotes the 
distance of the ion and the conducting wall in the case 
of the presence of a single wall. This result will need to 
be matched against a function D{Z), which measures the 
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effect in the case of two parallel slabs and sums over all 
mirror charges, with the proviso that the interaction here 
is proportional to the third inverse power of the distance. 
A series representation of the correction factor is easily 
obtained as follows, 


D{Z) 



1 

1 ) + Z]3 (an 



1 


1 



1 


12C(5) 

n3 


+ 0{Z‘^). 


2 

(an)3 
+ 4C(3) 
(43) 


The latter expansion confirms the consistency with the 
result for a single wall, as given in Eq. (42). The result 
given in Eq. (43) can be expressed in terms of the Hurwitz 
generalized zeta function. 


C(^a;) = 


n—0 


— C(3,x) = . (44) 

(n + x)^’ 2 


Eor perfect symmetry (ion perfectly aligned with the mid¬ 
point in between the two plates), the D coefficient is eval¬ 
uated as 


Dil,Z) = D(z=^) 


12C(3) 

Z3 


14.4247 

Z3 


< 2 X 


(45) 

The contribution of the additional mirror charges reduces 
the total result by about 10 % in comparison to the con¬ 
tribution from the two closest mirror charges alone, which 
would otherwise suggest a value of 16/a3 for the D coef¬ 
ficient. 


B. “Sniper” Setup Configuration 

The basic idea is obvious from Fig. 1. An ion enters 
the beamline, pre-accelerated. Energy selection (velocity 
selection) with a control measurement of the ion’s energy 
proceeds before entering the interaction region between 
the two aligned, parallel plates. For an accelerating volt¬ 
age of order 20 V, and single-charged helium ions (He"'"), 
the de Broglie wavelength is of the order of 5.1 x 10“33 ni. 
Thus, we can safely ignore diffraction effects which could 
otherwise occur when entering the interaction region in 
between the plates. Interactions with the image charges 
on both sides of the beam track can be ignored under per¬ 
fectly aligned (“sniper”) conditions. The idea is that the 
atom loses energy on its trajectory due to non-contact 
friction, with a corresponding measurable energy loss af¬ 
ter its has left the interaction region. An aperture before 
the ion’s energy measurement area, about a meter away, 
ensures that no significant distortion of the ion’s path 
due to interaction with the image charges has occurred. 
This is necessary because the Coulomb interaction with 
the image charges does work on the projectile ion, poten¬ 
tially altering its kinetic energy. 


The friction force in the interaction region is given as 

F\\ = -’7(2) Vx , (46) 

where r)( 2 ) is the damping coefficient for a configuration 
with two parallel plates, which we consider according to 
Eq. (42), within the approximation (45) for the ion in the 
middle of the parallel plates. The conductivities of met¬ 
als strongly depend on the concentration of dopants. We 
use the following values, which represent estimates of the 
room-temperature direct-current conductivity a = ctt(O) 
of a number of metals which can easily be formed into 
almost perfect slabs, either as a polycrystalline material 
with a well-polished surface, or even as single crystals (as 
in the case of vanadium), 

CTAu « 2.30 X 10^(Om)-i , crva - 5.08 x 10®(nm)-i , 

(47a) 

CTTi 1.27 X 10®(f2m)-i, crc 1.28 x 10^(Om)-i . 

(47b) 


Here, ac refers to graphite, where we assume that the 
basal plane of the hexagonal crystal lattice is aligned with 
the surface plane of the slabs. [The data in Eq. (47) has 
been compiled as the average of data given by manu¬ 
facturers for several commercially available metals with 
different polycrystalline structure and different dopants; 
in a precision experiment, it would seem indicated to 
measure the low-frequency (direct-current) limit of the 
conductivity cr of a specific sample independently.] For 
reference purposes, we here indicate that an evaluation 
of L according to Eqs. (29) and (47) for graphite leads 
to a value of L = 1.38 x 10“^® (rad/s)“^. For helium ions 
at a distance Z from the plates, with 

m = 6.646 X 10“^^ kg, Z = 0.5 ^m, (48) 


one has mdvx/dt = —’ 7 ( 2 ) Vx or 


dvx 

dt 


2’7(2) 

m 


Vx 


-Tv 


X ■ 


(49) 


Helium ions (He+) pre-accelerated to 20 eV energy en¬ 
ter the interaction region at a speed of 3.11 x 10^ m/s. 
Assuming a 10 cm long interaction region, the fractional 
loss in the flight velocity is calculated as 


r= l-exp(-rAt), At = 3.22 X 10"®s, (50) 


where At is the flight time in the interaction region. We 
obtain the following fractional kinetic energy losses, 

tau = 1-55 X 10“'^, rva = 7.03 x lO-"^, (51a) 

rxi = 2.82 X 10“® , re = 2.80 x 10■^ (51b) 


The corresponding relative decrease in the kinetic en¬ 
ergy (proportional to the square of the velocity) is twice 
as large as these values. The “sniper” aspect of the con¬ 
figuration sketched in Fig. 1 comes into play when we 
restrict, geometrically, the acceptance region for the ions 












leaving the interaction region. Let us assume that roughly 
Z? = 1 m further down the beam line, we restrict the 
available angular region for the arriving ions to a circu¬ 
lar aperture of diameter (see also Fig. 1) 

Os = 100 Z = 50 /rm. (52) 

For non-perfect alignment, the image charge interaction 
may exert a force in the interaction region, in the positive 
or negative z direction, without altering the x component 
Vx- The calculation of the interaction energy with the 
mirror charge(s) has been discussed in Sec. IV A. How¬ 
ever, in order to obtain a figure-of-merit for our proposed 
experimental setup, it is not sufficient to consider the ge¬ 
ometric requirement of passage through the geometric 
aperture indicated in Fig. 1; it corresponds to a restric¬ 
tion of the modulus of the velocity change according to 

where we redefined ^ as compared to Eq. (22). The figure- 
of-merit of the measurement is obtained as the ratio of 
the relative change in the velocity due to non-contact 
friction, divided by the geometrically restricted “uncer¬ 
tainty” in the velocity measurement, 

. 2(l-exp(-rAt)) 

C [as/DY 

We obtain 

/Au = 1.24xlO^ /va = 5.62 X 10^ , (55a) 

/Ti = 2.25 X 10^ , /c = 2.23 x 10^ . (55b) 

These results show that a measurement should be feasi¬ 
ble. We also have done an evaluation based on the exten¬ 
sive reference volume [23], analyzing available data for 
the dielectric function of a-quartz, along the ordinary 
(o) and extraordinary (e) axis. We obtain the following 
values for the L coefficient, 

Lo = 1.40 X 10“^’’ (rad/s)“^, (56a) 

Le = 2.14 X 10-1^ (rad/s)-i, (56b) 

which translates into the following loss coefficients and 
figures-of-merit. 

To = 2.82 X 10“® , Te = 4.30 X 10"® , (57a) 

/o = 9.04xl0^ /e = 1.37xl0^ (57b) 

A measurement using single-crystalline quartz surfaces 
also seems possible. The “sniper” configuration is “auto- 
correcting” in the sense that the ions can only pass 
through the detector aperture under well-aligned condi¬ 
tions. 


V. CONCLUSIONS 

We investigate the energy loss of an ion in the vicin¬ 
ity of two conducting surfaces, in a configuration where 
the electrostatic interaction with the two image charges 
compensate each other. Based on a rederivation of the 
frictional force, we relate the general expression for the 
friction force to the functional form of the permittivity at 
low frequencies, which we parameterize in terms of a sum 
of generalized Drude and Lorentz profiles. The definition 
of a conductor implies the existence of a zero-resonance- 
frequency term in the permittivity and we show that it 
is this term which dominates in the evaluation of the 
frictional force. 

Specifically, the result for the non-contact friction co¬ 
efficient, given in Eq. (21) for a single wall, is generalized 
to the interaction with two parallel walls in Eq. (42). 
The L coefficient can be written in terms of a parame¬ 
terization of the dielectric response functions of the ma¬ 
terial, according to Eq. (33). For a conductor, the limit 
L is exclusively determined by the conduction band [see 
Eq. (29)]. 

We identify the main problem in a conceivable exper¬ 
iment as the isolation of the frictional force from the 
strong electrostatic interaction with the image charge. 
In general, the functional form of the interaction poten¬ 
tials of ions with other electromagnetically interacting 
objects is different from the corresponding expressions 
for atoms. In Sec. I, we recall that the ion-atom interac¬ 
tion {1/RY is stronger than the atom-atom interaction 
( 1 /i?®), and that the friction force in the vicinity of a sur¬ 
face (1/Z^) is larger than that for an atom (1/Z^). In a 
typical case, there is a difference in the power law involv¬ 
ing two more inverse powers of the distance for “atom¬ 
something” interactions as compared to the correspond¬ 
ing “ion-something” interaction. 

In principle, the enhanced functional form of the non- 
contact friction coefficient for ions as compared to atoms 
would recommend a measurement of the quantum fric¬ 
tion effect using ions. However, the strong interaction 
with the image charge implies that a careful compensa¬ 
tion of the interactions with the images becomes neces¬ 
sary in the vicinity of a conducting surface, in order to 
minimize a deflection of the ion’s trajectory. We aim to 
identify a collection of suitable parameters, in terms of 
a pre-accelerating voltage, to find a compromise between 
a slow flight velocity (which enhances the “risk of de¬ 
flection” ) and a too short interaction time (which would 
make the quantum friction effect undetectable). A pro¬ 
posed set-up with certain auto-compensating features, to 
isolate the friction effect, is described in Sec. IV. The 
most suitable materials for the experiment are be those 
which are readily available as large-scale single crystals, 
and constitute conductors, but with a low value of the 
conductivity, in order to increase the Ohmic heating due 
to the flow of the image charge. We find that the best 
figure-of-merit is obtained for graphite, while a readily 
available metal like vanadium, which has excellent sur- 
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face properties, also would be available for a precision 
experiment. 
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Appendix A: Fourier Transform in Signal Theory 

In signal theory, different conventions are used for the 
Fourier transform. Normally, in physics, one integrates 


the Fourier frequency transform with the “integration 
measure” dw/(27r). However, according to the Eq. (3.118) 
of Ref. [3], and Eqs. (1.90) and (1.91) of Ref. [4], we have 

f(t) = J dwe”'"^*/(w), (Ala) 

/M= / ^e-‘/(t), (Alb) 

and the same conventions are employed in signal the¬ 
ory for the corresponding transformations from coordi¬ 
nate space to wave number space, namely, one integrates 
with the integration measures d^r/(27r)^ and d^k instead 
of d^r and d^fc/(27r)^. The difference in the integration 
measures also explains the occurrence of multiplicative 
factors of 27r in other formulations of the fluctuation- 
dissipation theorem, such as Eq. (3.118) of Ref. [3]. 
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